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Abstract 

We consider the contribution to the three-point function of matter density fluctuations from nonlinear 
growth after modes re-enter the horizon, and discuss effects that must be included in order to 
predict the three-point function with an accuracy comparable to primordial nongaussianities with 
/nl ~ few. In particular, we note that the shortest wavelength modes measured in galaxy surveys 
entered the horizon during the radiation era, and, as a result, the radiation era modifies their three- 
point function by a magnitude equivalent to /nl ~ C(4). On longer wavelengths, where the radiation 
era is negligible, we find that the corrections to the nonlinear growth from relativistic effects become 
important at the level /nl ~ few. We implement a simple method for numerically calculating the 
three-point function, by solving the second-order equations of motion for the perturbations with the 
first order perturbations providing a source. 



1 Introduction 

The spectrum of initial density perturbations from inflation is known to be nearly gaussian. 
Should deviations from a gaussian distribution be detected, however, they will imply poten- 
tially powerful constraints on models of inflation. In particular, the shape and size of the 
three-point function of primordial density perturbations in many models is predicted to be 
near current experimental limits, and is one of the best observables for distinguishing different 
inflationary models pj^[3l|4j[5],[6l[7j[8],[9l EEl [IU EG] as well as its (perhaps less compelling) 
alternatives [13J. The best measurements so far come from cosmic microwave background 
(CMB) data [HI [151 HE], but large-scale structure (LSS) measurements are improving and 
will become comparable. The shortest wavelength measured by LSS is smaller than for the 
CMB, and it is a three-dimensional map instead of a two-dimensional one. Therefore LSS 
measurements potentially include many more modes. However, they also therefore involve 
modes that entered the horizon at earlier times, when the universe was radiation dominated. 
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Our purpose is to quantify the effect of radiation on nongaussianities in LSS measurements 
in comparison to that from primordial nongaussianities. 

The significance of the radiation era is due to the nonlinear growth of nongaussianities, 
which occurs even for gaussian initial conditions [T71 [18]. The leading contributions have 
been studied in the perturbation theory (PT) formalism, dropping the radiation component 
of the universe and using Newtonian gravity [191 |20~[ |2T] . While this approach has been very 
successful, the future improvement of observational data will necessitate knowing the theoret- 
ical error arising from the PT formalism assumptions. In fact, the late-time contribution to 
nongaussianities from PT is orders of magnitude larger than that expected from primordial 
nongaussianity [22], and even relatively small corrections to the PT ansatz may swamp the 
signal we are interested in. 

Experimental limits on non-Gaussianities are generically given in terms of a scalar variable 
/nl [2"3"ll2~i] which in the so-called "local" ansatz parameterizes the deviation of the Newtonian 
potential from a gaussian variable $ 5 as 

There are other possible shapes from models of inflation, but we will use this one for compari- 
son. We find that radiation effects on the LSS are comparable in size to /nl ~ 4. Current LSS 
experimental limits are still consistent with a gaussian distribution and have 95% confidence 
limits of —29 < /nl < 69 [25]. So far, they are nicely consistent with the limits from the CMB 
(26], which, when combined together, give the constraint —1 < /nl < 63 at 95% confidence 
level. 

2 Estimates 
2.1 PT Formalism 

We are ultimately interested in calculating the difference between the three-point function 
from the PT formalism and from a universe with a radiation component, and comparing this 
difference to the three-point function from primordial nongaussianities. Consider first what 
kind of effects we expect for small vs. large wavelength modes. For very large scales that cross 
the horizon deep in the matter dominated region, we do not expect any growth in 5 before 
the matter-dominated region begins and thus we can essentially treat the entire history of 
the universe as matter dominated. However, for short wavelength modes, growth has already 
begun deep in the radiation dominated region and it is not immediately clear whether this 
can give a large correction or not to the distribution of perturbations. The perturbations 
are very small during the radiation dominated era, and this implies that non-gaussianities 
generated during this region will be very small. That is, non-gaussianities from gravitational 
instability are generated because the first order fluctuations are a source for the second order 
fluctuations: 

5® ~ (2) 



2 



and thus this generates non-gaussianity dominantly at recent times. In this section, we will 
estimate the effect of radiation by considering how a large distortion in 5 during the radiation 
era affects the three-point function at late times. 

We first recall a few results. The method of PT is essentially to perform an expansion in 
the scale factor a. At second order, 

5(a,k) = 5i(k)a+ \5 2 (k)a 2 , (3) 
d 3 q 1 d 3 q 2 



82(h) = j d ^^ 2 5(k-q l - q 2 )F 2 (q 1 ,q 2 )5 1 (qi)5 1 (q 2 ) ■ 
In Newtonian gravity with matter only and no decaying mode, one obtains 



(4) 



\F 2 { qi ,q 2 ) = U^(^ + q ^)+^q l -q 2 f. (5) 



7 2 \q 2 q 1 J 7 

Thus, as a function of a, the three-point function is (5^5^ 5^) (a) ~ a 4 P 2 F 2 + a 3 Bi(ki, k 2 , k 3 ), 
where -P(-Bi) is the two-point (three-point) function of 5i, □ and thus Bi encodes the three- 
point function initial condition. Modifying the evolution history up to a scale factor around a eq 
will modify B\ proportionally by a cq . For instance, suppressing all growth of nongaussianity 
in this simple approximation up until a eq , so that (5^5^^) (a eq ) = 0, imposes 

B 1 ~ -a eq P 2 F 2 . (6) 

We therefore ought to expect corrections to the three-point function today proportional to 
o eq . 

Such a small correction matters because at late times and short wavelengths, the nongaus- 
sianity from nonlinear growth (i.e. from F 2 ) swamps the contribution from initial conditions. 
Within linear theory, 5k (a) = Mfc(a)$£ nm , 

M k (a) = - 3 p^§ Dl (a), (7) 

where D\{a) is the growing mode for 5 and T(k) is the transfer function. Under the local 
ansatz $P nm = $ 9 + /nl^, this leads to a dark matter three-point function proportional to 
/nl : 

B L {k X: k 2 M) = M fcl M fe2 M fc3 2/ NL (P $ (A; 1 )P $ (A ;2 ) + cyc.) . (8) 

The full three-point function is then approximately the sum of the contribution Bi from 
linear growth and the contribution Bq from non-linear growth. For the sake of this ap- 
proximation, we will take isoceles triangle configurations \k\\ = \k 2 \ = y\k^\, in which case 



1 The reader familiar with the PT literature should note the additional factor of \ in the expression for Fi. 
This is a consequence of the additional factor of ^ in our convention for 5^ 2 \ 

2 Power spectra for fluctuations other than S will be denoted with a subscript to indicate the fluctuation, 
e.g. P^{k) for <f> fluctuations. 



3 



\F 2 (k 1} k 2 ) = 10- and \F 2 {k x ,h) = ±F 2 (k 2 ,k 3 ) = £(13 - £). Note that y > \. We 
will also ignore the scalar tilt n s , so P$(k) oc k~ 3 . Now compare this to a correction in the 
nonlinear bispectrum of the size we expect from modifying the model at scale factor a eq , i.e. 
SB G ~ a eq F 2 (A;i, k 2 )P(h)P(k 2 ) + cycl.: 



/ 

SBq a eq M k3 
Bl ~ 2f : 



NL 



2y 3 F 2 (k 1 ,k 3 ) + F 2 (k 1 ,k 2 ) 
2y 3 + 1 

V / 



(9) 



We can simplify this in various limits. The transfer functions T(k) simplify at large and at 
small k in a flat (matter + radiation) universe: 

T(k « A; eq ) = 1 , (10) 

> 15A; cq ) = l g(A;/6A; cq ) . (11) 

z /c a eq 

In the limit of squeezed triangles y ^> 1, we arrive at the following estimate 

SBo w ^13 = f ^ (i) 2 y>>1 
Bl /nl 28 | ^y]og(*s/6*J fc 3 » fceq ' 

Note that all dependence on the shorter wavelength modes k±, k 2 drops out. Thus, even if k\, k 2 
enter the horizon deep in the radiation era, the 3-pt function is unaffected by the radiation 
corrections as long as k 3 enters during matter dominance (MD), where SBg is suppressed by 
(k 3 /k eq ) 2 ; however, the 3-pt function is affected if k 3 enters during the radiation era. 
We next turn to the limit of equilateral triangles, y — l,k = kf 

**» „ t^Mil = { sS£ (£)' k « k «> i y = i (is) 

Bl /nl 7 1 =M hg(k/6K,) 

In other words, modifying the matter content of the universe at a eq near or after matter- 
radiation equality alters the bispectrum roughly corresponding to /nl ~ few for modes that 
entered the horizon somewhat before matter-radiation equality. We will need a more pre- 
cise calculation of nonlinear growth near matter-domination equality in order to predict the 
bispectrum at order / NL ~ few. 



2.2 Radiation Effect in Newtonian Gravity 

By the previous arguments, we expect that the dominant contribution from radiation will 
come from the era near the transition between matter and radiation, because at this time, S 
is as large as possible before matter dominates. Short wavelength modes, which by previous 
arguments will get the largest contribution, will therefore already be inside the horizon at 
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that point and we should be able to approximate the effect on them using Newtonian gravity. 
In this limit, we will be able to obtain an analytic approximation. 

The equations of motion for the dark matter perturbations in Newtonian gravity are (e.g. 

CHI) 

S' k + ikV k = - (l + ^\ S qi iq 2 V q2 , (14) 

ikVi + HikV k + ^H 2 n m (v)$k = \k 2 qx ■ q 2 V qi V q2 . (15) 

Primes denote derivatives with respect to 77. The behavior of the background is that of a 
universe with matter and radiation with no cosmological constant, which can be written in 
terms of a eq = ^f- as [27] 

a o 1 2 2 H o V2-1 . . 

— = 2ar] + a r\ , a = — = , (lb) 

°eq 2 A /a eq (l + a eq ) Veq 

and VL m {rf) = = -n . It is straightforward to solve equations (114|15p at linear order. 
The solution [j normalized so that 5^ = 1 today is 

ff)fo )= 2 + 3^(2 + ^) lW (i) () = Kl+^j 

9 v/; 2 + 3«r/ (2 + ar] ) 9 yu 2 + 3ar/ (2 + ar) ) K ! 

Taking the first order solutions to contain only a growing mode, we obtain an analytic solution 
for the second order modes, which for brevity we do not write out. Instead we expand in 1/a, 
and obtain the PT result equation ((Sj) plus corrections. The leading correction is enhanced 
by log(a eq ): 

\F 2 ,rad = \F 2 , PT - |a eq log (a^ 1 ) \ 2 (k 2 - (qi + q 2 f) (k 2 - (q x - q 2 ) 2 ) + ... 

4 

= \F 2 , PT + |a eq log (a" 1 ) — (1 - (ft ■ q 2 ) 2 ) +... . 



For example, on equilateral configurations (i.e. k — qi — q 2 ), we have F 2ira d — F 2y PT ~ 
a eq \og(a~q)F 2PT , which is ~ 1.2a eq F 2PT for a eq = 3 x 10~ 4 . As in the case without 
radiation, the kernel F 2;T . ad is invariant under rescalings (k, q±, q 2 ) — > X(k, qi, q 2 ), and all scale- 
dependence of 5^ comes from the scale-dependence of 5^. Note that the leading term 
vanishes on "squashed" triangle configurations k — \q\ ± q 2 \, and in this case the sub leading, 
non-log-enhanced contribution dominates. However, the non-log-enhanced contribution is 
more sensitive to the behavior of the linear solutions during the radiation era, and in order 
to get an accurate estimate valid at large k we would have to include some effects we have 
neglected. 



3 We take the initial conditions at arbitrarily small 77, so the decaying mode is projected out and the solution 
is the growing mode. 
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In the above approximation, we included radiation only in the background dependence, 
i.e. in TC(t]). This approximation is exact in the limit k 3> k eq . A more careful treatment 
would include the effect of the radiation density perturbations as well, which would enhance 
the evolution of the Newtonian potential $ at early times. One might then obtain an analytic 
approximation by including the backreaction of only the radiation on the Newtonian potential 
at early times and only the matter at later times, and matching the two solutions in an 
intermediate regime. One would expect in this way to see an enhancement of the final 
perturbation size, similar to results from the analogous method in linear perturbation theory 
|28j . However, we instead will now turn to a numeric calculation that will give the fluctuations 
at second order for all wavelengths. 

3 Method 

In this section we will describe our method for calculating the effects of radiation on the dark 
matter density three-point function at late times, based on an expansion at second order in 
the fluctuations. That is, we split all perturbations O = + into a first order piece 

and a second order piece 0^ 2 \ For closely related work on second order corrections 
to cosmological perturbations see e.g. [29j USB EB E21 [33j [Mj. Since we are interested in 
following the dark matter distribution, we can approximate the baryon-photon fluid as a fluid 
with w = |, and otherwise neglect Qb for simplicity, taking a universe of two perfect fluids 
composed respectively of collisionless matter and radiation. In this approach, the second 
order perturbations are sourced due to nonlinear interactions by the first order ones, and thus 
are quadratic in them, ~ O^O^. This leads to a nonvanishing three-point function 
even for purely gaussian initial conditions, and this contribution should be removed in order 
to reconstruct the primordial nongaussianity. In particular, we are interested in the kernel Fi 
in Fourier space: 

¥f^) = \ J -qi- &)F 2 {qi, <f 2 ; v)5&(vo)s&{no) , (19) 

where rj is conformal time and we choose i]q = rftoday in order to make contact with existing 
literature. We compute F 2 in a straightforward manner for any qi,q 2 by inputting the first 
order solutions into the equations of motion as sources and solving numerically for the second 
order terms. This procedure could be carried out iteratively to higher orders as well by 
inputting the first and second order solutions as sources at third order, etc. 

We also may neglect all vector and tensor degrees of freedom. The reason is that vector 
and tensor degrees of freedom vanish at first order in the initial scalar perturbations, arising 
only at second order. Therefore, we may decompose any equation with a vector index into its 
scalar O s and vector % v piece. In Fourier space, the most general equation at second order is 

{aq\ + bq\)0^\ qi )OV{q 2 ) + VOf\k) + Of\k) = 0, (20) 
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where k l O l v (k) = 0. By contracting this equation with k l , we therefore project out all vector 
modes and obtain a second order equation for just the scalar modes. A similar argument 
applies to the tensor fluctuations that arise from primordial scalar fluctuations. In addition, 
tensor modes have their own primordial fluctuations, whose size is typically suppressed with 
respect to the scalar fluctuations in models of inflation. More important, their two-point 
function with the scalar modes necessarily vanishes by rotational invariance, and thus they 
give no contribution to the matter three-point function. We may therefore neglect them as 
well. 

The computation is simplest in conformal Newtonian gauge, where the metric is diagonal. 
Specifically, we will take 

ds 2 = a 2 (-e^dT] 2 + 2u i dx i dr] + (e" 2 % + X ij)dx l dx j ) , (21) 
u*i = 0, XW = 0, (22) 

so that uJiiXij contain no scalar modes. We may therefore neglect u and \ completely for 
following the scalar modes at second order. 

The energy- momentum tensor is that for a fluid of matter, which is pressureless up to 
negligible O(-) corrections, and a fluid of radiation w = 1/3: 

(A 1 \ /e - * v 2 mT e* i 

T^ u ,m = Pm U m U v,m j ^t,r = Pr ( 7^ U r U v,r + J > U m,r = I 0- ~l ^ )) ~^ Vm ^ r 

(23) 

We will denote the longitudinal piece of velocity as v l = k l V . We obtain the equations of 
motions from Einstein's equation G lL v = 87rG N T^ and the conservation equations V^T^ mr = 
0, V ^(Pm^O = 0, which is valid because matter and radiation interact only gravitationally. 

It is straightforward to derive the Einstein tensor exactly in conformal Newtonian gauge 
with only scalar modes. We write it down for reference in appendix [A] There are two 
components of the Einstein equation that will be useful to us. The first is the time-time 
component, a 2 G° = 3H 2 T° Q /p, which in Fourier space is 

k 2 ^ 2) + 3w*J a ) + m 2 ^ + \n 2 ( ^4 2) + ^ 2) ) = -5 3 = (24) 

2 V P P J 

(?i ■ Q2 - 4g 2 2 + m 2 ) + + 

— 37i 2 qi ■ (J2 ( ^-V qi m V q2m + ~~Vqi ,rVq 2 ,r 

VP 3 p 

Primes denote t] derivatives, and TC = In the above equation, an integral (2tt)~ 3 J d 3 qid 3 q 2 5(k— 
Qi ~ 92 ) is implied over terms quadratic in the first order perturbations. In general, we mean 
for such an integral to be implicit in any equation in Fourier space with terms quadratic in 
the first order perturbations. The second useful component is the shear piece, specifically 
{hy — ^5ij)G l j, which projects out the 5^ piece in £?* .. The resulting equation of motion is 

1 -k 2 (-^ + ¥V) = 1 -k 2 S 4 (k 1 r ] )= (25) 
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((-2(k. qi )(k.q 2 ) + ~ qi . q2 )$U<!>U 

+m\(k ■ q x ){k ■ q 2 ) - ift • q 2 ) (^V qi>m V q2iTn + ^V qi , r V q2 , r 

We now turn to the conservation equations for matter. The first comes from V /J (p m M^ l ) = 0, 
the second from V^T^ m = 0: 

SS i +ikVS ) -3^=S 1 (k iV ) = 

^^W^+^v|,( 2 ) = S 2 (M) = 

2(* • &)$f'y£L - • ^X.l^L - 2i*$g)*g) . 

where we have simplified the sources Si by imposing the first-order equations of motion and 
symmetrizing in qi,q 2 when convenient. 

We take the two radiation equations of motion from V^T^ r = and V M T^ r = 0: 

ifO _ 4$< 2 ) + = S 6 = (28) 

g •"dl Hi Y qi,r 12,T g ^ x " 2 gi,r 92 g " z qi Q2,r ' qi,r (}2,r ' 

K' (2) + j^ 2) + ^ (2) = 5 7 = (29) 



^ {2q iq2 + {q\ + q\ - 3k 2 )q, • &) V#> V® + £*£^£>, + 4g 2 • IfW^ 



The behavior of the background is that of a universe with matter and radiation, as in 
equation ffTET) . In addition to the above equations of motion, we require the equations for 
the first order perturbations. The equations of motion for the first order perturbations are 
identical in form to those for the second order solutions, with the sources Si set to zero. We 
solve these numerically. 

The initial conditions for the fluctuations depends on physics before they reenter the 
horizon. There is both a contribution from before a mode exits the horizon during inflation 
that depends on the inflationary model as well as a possible contribution due to light degrees 
of freedom that perturb the reheating surface. Outside the horizon, it is useful to work with 
the metric of scalar modes can be set by comparing with ( gauge, which in the absence of 
spatial gradients is defined at nonlinear order by 

ds 2 = a 2 (rj) (-drj 2 + e 2C dx i dx i ) . (30) 

For the minimal model of inflation, ( is initially an approximately gaussian field with a three- 
point function that is suppressed by a slow-roll parameter [TJ [21 [35]. We therefore choose 
£( 2 ) = initially, in order to separate out the effects due to nonlinear growth at late times. 
This condition in turn sets the initial conditions on most of the remaining perturbations 
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through the equations of motion. Let us first consider how it sets and $^ 2 ^ initially. 
and are related to each other through the constraint equation (I25p . which enforces 

^(2) = $ (2) + ^ _ (31) 

At nonlinear order, transforming from ( gauge (I30p to Newtonian gauge shows that ( = 
— $ — during the radiation era in the absence of spatial gradients. Thus, at a -C a eq , kr\ <C 1 

* = ~C - \Sa , (32) 

6 

and therefore, we have |$^ 2 ^ = — -54 initially. 

The initial conditions for the perturbations therefore depend on the behavior of the source 
terms at early times a <C a eq during the radiation era and long wavelengths kr\ -C 1, because 
as above the equations of motion each contain perturbations without time derivatives. At 
early times, these perturbations contain only a growing mode, which is usually constant 
outside the horizon, and the equations of motion turn into an algebraic relations between the 
perturbations and the sources. This is also what happens in the linear theory, except that the 
source terms Si are not present in that case. In the long- wavelength early-time limit, we will 
see that we need to use only S 3 — > — QTi 2 ^^^ and S4 — > — p- (k ■ q\k ■ q 2 — \q\ ■ q^j 

(using iVr 1 ^ = |A;77$^ initially). Equation (12^|) and ([25]) together imply the initial condition 
for 5i 2) : 

S ® = ~^ S3 " " 2$(2) ' (33) 

(2) (2) 

The second order velocities Vm , V r , start off negligible at early times; if needed, they can 

be read off of equations f)27|29|) by taking oc i] at small 77. 

(2) 

Finally, we need to set the initial condition for 5m ■ It does not appear in the equations of 
motion at early times, because it always has a time derivative acting on it or else is multiplied 
by (p m /p), which vanishes in the infinite past. We will assume adiabatic initial conditions, 
where there is just a single light scalar degree of freedom outside the horizon (i.e. no entropic 
modes). In this case, we may set the initial condition on 5$ by the fact that all scalar 
fluctuations arise from the time-shift of their background. More precisely, for adiabatic initial 
conditions, all scalar fluctuations outside the horizon arise from a single scalar fluctuation, 
which may be parameterized as the pion 5t for the spontaneously broken time translations 
[Sj. Thus, for either matter or radiation, p + Sp = p(t + St) = p{i) + 5tp(t) + ^5t 2 p(t) + . . .. At 
first order, 5t^ = ^ is therefore equal for matter and radiation, and this implies the usual 

adiabatic relation 5 m = |5 r initially. At second order, one must equate 5t^ = — | (~^~) > 
and therefore the initial condition for dm is 
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In principle, the equations of motion (I24H29P and the initial conditions are all that is needed 
to find the second order kernel on any scale. In practice, however, we also must specify 

(2) 

the gauge in which o m is actually measured in observations. While it is straightforward 

(2) 

now to calculate 5 m in conformal Newtonian gauge, it is unlikely that this is the quantity 
we are interested in. We shall not determine in this paper which is the correct gauge for 
5 m corresponding to measurements (or, put differently, what is the correct gauge-invariant 
observable measured in surveys). Rather, we shall treat gauge dependence as an additional 
source of uncertainty. What this means requires some care, as one can always choose the 
gauge 5 m = 0, where the matter density provides the clock for the coordinate system. We 
shall avoid such extreme gauge choices and instead consider the change in 6^ going from 
Newtonian to synchronous gauge. On subhorizon scales k/Ti ^> 1, perturbations become 
increasingly local and insensitive to the gauge choice, and S from one gauge to the next 
changes as ~ ^5. For example, under a time diffeomorphism r\ — > r\ + a, \1/ changes at first 
order by ^ — ► ^ — Tia — a', and 5 changes at second order by 5 — > 5 + 3Ha(l + 5) + . . .. Thus 
inside the horizon during the matter era, a gauge transformation that changes the metric \I> 
by 0(1) has Ha ~ and consequently 5< 2 > ~ 3^5^ ~ 3^4^ ~ ^{^f. Thus the 
physical significance of contributions to F 2 that behave parametrically as 7i 2 /k 2 is unclear, 
and much additional work is required to understand these contributions. 

4 Results 

4.1 Numeric Results 

The second order kernel F 2 is a function of two three-momenta <fi and q 2 and is thus naively 
a function of six variables. However, due to rotational invariance, it is in fact a function of 
only three variables, which we will choose to be k, X\ = y, and x 2 = In the limit of short 
wavelengths and negligible radiation, F 2 approaches the PT result of equation (jSJ), which is 
scale invariant and thus depends only on x\ and x 2 . Note that 1 + x 2 > X\ > 1 — x 2 by 
momentum conservation and we may take x\ > x 2 without loss of generality by symmetry. 
The kernel F 2 (q\,q 2 ) enters the three-point function through 



and in the three-point function one may further take ki > k 2 > k^ without loss of generality. 
It is convenient to introduce the reduced three-point function 



(6 kl S k Jk 3 ) = S [J2 ki) (2tt) 3 (F 2 (h, k 2 )P kl P k2 + eye.) 




(35) 



Q{k 1 ,k 2 ,k 3 ) 




(36) 



P(h)P(k 2 ) + P(k 2 )P(k 3 ) + P(h)P(k 3 ) ' 
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Numeric 



1.0 




5g(xl, x2) 



Figure 1: The kernel x l x 2 1 Q(x\, x^) for k fixed at k 



11/ceq- The left plot is from our 



numeric computation, and is indistinguishable by eye from the PT limit given in equation 
(J5j). The right plot is the difference between the numeric result and the PT result, rescaled 
by a eq , because by equation (jSD, a eq controls the size of the radiation correction. 

Q(ki, k-ii k-&) is identical to F 2 (k, k) on equilateral triangles. Unlike F 2 , though, it is symmetric 
under permutations of ki,k 2 ,k%, and thus takes into account some cancellations in the three- 
point function between F 2 (ki, k 2 ) and its permutations. 

A plot of x\ x x 2 Q(x±, X2) from our numeric computation with k fixed at 11 k eq is plotted in 
Figure [IJ along with a plot of the difference between our numeric result and the PT result @. 
The observed value of k eq is 0.013/iMpc _1 . At k = 11 k eq , the difference is very close to 
that obtained in Newtonian gravity with radiation, equation (1181) . shown in Figure [2J This 
indicates that > 10 is already large enough that the approximation used in section 12.21 is 



The result from our numeric method takes into account not only radiation corrections 
but also corrections from general relativity (GR), the leading contributions of which are post- 
Newtonian (PN) corrections. The scale-dependence of these two effects is completely different. 
GR effects are important on larger scales and behave parametrically as (7i/k) 2 compared to 
the PT kernel. We can partially isolate such effects by using our numeric method with 
Q r = (equivalently, a eq = 0), and taking the initial conditions appropriate for a matter-only 
universe. In fact, in this case the second order perturbations may be solved for analytically, 
as in [301 ED El], an d seen explicitly to be suppressed by (H/k) 2 with respect to the leading 

4 We define the "difference" 8Q as the change from corrections to F2, and not from corrections to the 
power spectrum P(k). Explicitly, 



good. 



^2(91,52) 



SQ 




(37) 
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Analytic Difference 




Figure 2: The log-enhanced difference between the PT kernel and the kernel from using 
Newtonian gravity with a background radiation component. The analytic result is given in 
equation ([TBI . 



piece, as shown in appendix lAl (eq. l60l) : 



2°k 



(Pit - Ofc) + • ?2j I — + — I + a&(gi • g 2 ) + 7fc 

2 \q 2 QiJ \Q2 Q2 



S 0) gW 
u qi u q2 

?r it 



2 59ft 2 45ft 4 _ ft 2 54ft 4 _3ft 2 9ft 4 

~ 7 + MP + ~2F~ ' ^ ~ 2fc2 + ' 7fc ~ ~ 2P + 2F " 

As we have discussed, however, such effects are gauge-dependent and require more care 
in order to be written in terms of present-day observables. The effect of radiation, on the 
other hand, is nearly independent of the magnitude of k for modes that entered the horizon 
in the radiation era and approaches a constant correction at small scales. In order to give 
a sense of the scale-dependence, we may fix the shape of the triangle formed by k, qi, and 
§2- In Figure [3], we take ki = k% = yk, k\ = k and plot the correction to Q as a function of 
k/H. For comparison, we also show the correction in a matter-only universe as well as the 
leading correction (Tl8|) in Newtonian gravity when radiation is included in the background 
(but radiation perturbations are neglected). For our choice of parameters, k eq /H m 80, and 
the radiation correction starts to dominate over PN effects between about k eq and 10 k eq , 
depending on the triangle shape. The small parameter suppressing the PN effects is (H/k) 2 , 
and thus they are larger than the radiation effects when (H/k) 2 > a eq ~ (H/k eq ) 2 , i.e. when 
k < A; oq . The radiation effects and PN effects have opposite sign for equilateral and squeezed 
triangles, and thus the total correction crosses through zero in the transition region. One 
can see that the second order perturbation 5^ approaches our simple approximation (eq. 
[THjl . In fact, from Figure HI we see that the corrections over the entire range of k are well- 
approximated by the sum of the radiation corrections (eq. [T8|) and the GR corrections (eq. 
160]) : 5Q2 ~ 8Q2,rad + 8Q2,gr- Note that, although the radiation correction to Q2 grows as a 
function of redshift a little slower than (1 + z), the correction from primordial nongaussianities 
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y=l/2 y=l y=10 




k/H k/H k/H 



Figure 3: The correction \5Q\ for k\ — k^ — yk^^k^, = k with y fixed at |, 1, 10, from left 
to right, respectively. The correction is shown for a) a matter-only universe (a eq = 0), red, 
dotted, b) the exact numeric correction, black, solid, and c) the difference in Newtonian 
gravity between a universe with matter and one with matter + radiation (blue, dot-dashed, 
for which the log-enhanced piece of equation ([18]) vanishes when y — |; the size is independent 
of k because it is the difference between two results within Newtonian gravity). The black 
dashed line is the exact radiation effect, defined as the difference between the exact numeric 
correction and the matter-only correction. The exact numeric calculation has a eq = 3 x 10~ 4 , 
in which case k eq /H rs 80. 

has nearly the same z dependence and their ratio is nearly flat, as shown in Figure [5] 

Finally, we show the comparison between the corrections to 5Q and the contribution 
from primordial nongaussianities with /nl = 1 in Figure [6j The correction from radiation 
becomes the dominant correction only around k > 3fc eq - At k = 10k eq = 0.13/iMpc -1 , which 
is comparable to the largest k that next-generation LSS observations will be able to resolve 
[TTl [22] . the correction from radiation is comparable to /nl ~ 4.6 for equilateral triangles. 

4.2 Including A 

So far, we have considered a universe with Q m + Q r = 1, which is only a good approximation 
up until grows to be non-negligible. To calculate the effect of radiation in our universe 
with vacuum energy A ^ 0, one simply includes the effect of A on the background, which 
modifies 7i (77). Otherwise, the equations of motion for the linear and second order solutions 
are unchanged. 

In a universe with matter and A, the Newtonian equations of motion do not admit an 
exact solution which is separable, i.e. of the form 5^ (77) ~ F 2 (qi, <?2)^i(<?i)^i(?2)(-D+^( r ?)) 2 - 
It has been argued that, to good approximation, they are almost separable (see e.g. [18J), 
and D+(j)) takes the form of the (matter + A) linear solutions whereas -^(^i, Q2) is the same 
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Figure 4: Comparison between exact numeric corrections and an analytic approximation 
8Q 2 = 8Q2,rad + $Q2,gr- The corrections \SQ\ for k% = k 2 = yk 3 = yk with y fixed at ~, 1, 10, 
for left, center, right, respectively. In the top row, the exact result (analytic approximation) 
is shown in solid, black (blue, dashed). The difference between the exact and approximate 
results for SQ are shown in the second row, divided by their sum in quadrature. 



as the PT result. A better approximation is to calculate 6 k exactly in Newtonian gravity, 
matching the linear solutions onto the matter era growing modes 5 oc a at a <C 1. When we 
compute 5^ numerically in this way, we find that the result is exactly matched by the form 

lF 2 ( gi , q 2 ; Q m ) = e(Q m ) + ^ + e (n m )) (& ■ q 2 f , (39) 

2 V<?2 qij 

where e(Q m ) is determined numerically. A parameterization of e(fl m ) that turns out to work 
well is 

e ex act(On) « e flt (O m ) = ^ + c(l - fi^ /15 ) , c = 0.023969 . (40) 

As we show in Figure [7J this parameterization fits the exact numeric e(Q m ) to better than 
0.0035% for 0.2 < Q m < 1. 

The exact Newtonian result provides a much better approximation to the second order 
density perturbations than the PT result does. Since the linear modes are normalized at 
their current values, the PT three-point function today is insensitive to the value of fl m . For 
comparison, in Figured] we show the difference between the PT result and the exact result for 
the reduced three-point function, as well as the difference between the Newtonian result and 
the exact result. Including A, the Newtonian result is still only off by about the equivalent of 
/nl ~ 4, due to radiation. However, the unmodified PT result of equation ([5]) is off by over 
an order of magnitude more. 
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Redshift z Redshift z 



Figure 5: The correction \SQ\ relative to the primordial contribution with /nl = 1 on equi- 
lateral triangles as a function of redshift z. The wavenumber in the left plot is, from top to 
bottom, k = k eq x {30, 10, 4, 3, 2.3}. As k decreases from 30fc eq , the GR corrections grow and 
approximately cancel the radiation corrections near k ~ 2k cq . At smaller k, the correction 
from GR is well-approximated by the matter-only analytic result fl59l) . The right plot shows 
the exact result relative to 5Qf NL= i for k = 0.4/c eq , where essentially all the correction is from 
GR. The reason the correction in the left plot decreases at large z is that 5/<p gets smaller in 
the past, so the contribution from /nl is bigger in terms of 5. 




Figure 6: The corrections \5Q\ for k\ = k 2 = yk 3 = yk with y fixed at |, 1, 10, for left, 
middle right, respectively, compared to the contribution from primordial nongaussinities with 
/nl — 1- The exact numeric correction to the PT kernel is dashed, whereas the contribution 
from /nl is shown in solid. 



4.3 Leading Effect of Non-Gaussianities on Galaxy Power Spec- 
trum 

Ultimately, the dark matter three-point function must be related to observables. In [22], it 
was proposed that non-gaussianities in the dark matter distribution might best be measured 
through the additional scale-dependence they induced in the galaxy two-point function. This 
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Figure 7: Comparison of the exact numeric result for e(O m ) compared with the parameteri- 
zaton 0.023969(1 - fi™ /15 ) for a fiat universe with A and matter only. 
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Figure 8: The corrections in a flat universe with fi m = 0.27, f2 ri0 = 3 x 10.~ 4 Q m to |&3| 
for k\ = k 2 — yk 3 = yk with y fixed at ~, 1, 10, for left, middle, right, respectively, compared 
to the contribution from primordial nongaussinities with /nl = 1- The black, solid line is 
the difference between the exact result and the PT result. The blue, dot-dashed line is the 
difference between the exact result and the Newtonian approximation of equation fl39l) . The 
contribution from /nl = 1 is shown in dashed, black. The straight red, dotted line is the 
Newtonian radiation correction of section 12.21 including the non-log-enhanced piece. The PT 
result is very good on squashed triangles, but off by more than an order of magnitude more 
than the Newtonian approximation on equilateral and squeezed triangles at large k. 



approach has the advantage that it is possible to write the effect in terms of quantities that 
may be studied with exactly gaussian statistics. 

The reason for the leading effect on the two-point function is that non-gaussianities cor- 
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relate the long- and short-wavelength fluctuations. With gaussian fields, a long-wavelength 
fluctuations looks like a background fluctuation that increases or decreases the amount of 
short- wavelength fluctuation needed to form a bound object [3H]- However, when differ- 
ent modes are correlated, a long-wavelength fluctuation tends to increase or decrease the 
actual size of the short-wavelength fluctuation. More precisely, with <C ks, 8 kh is a long- 
wavelength mode that affects the variance of the short- wavelength mode 5k s : 



(4,4,) = ((5F 2 (h, k L ) + -rj^)P fel + (^ 2 (fc 2 , fci) + tt^)^) x5(h + k 2 + k L ) 

-k L Mk L Mji L V 



(41) 



where the subscript means that the expectation value is taken with 8 kL equal to zero, and 
where we have included the contribution from / NL through 8$ = M ks (^^ +f N j^(^^~k^^) ks ) 
in linear perturbation theory, with the * indicating a convolution. We are also concentrating 
on the difference between the effect due to radiation and general relativity with respect to the 
standard newtonian treatment. The value of the short-wavelength mode variance determines 
the local value of the rms density fluctuations, i.e. og, and thus the total effect of the long 
wavelength perturbation is (see for a more precise derivation [25, 39J): 

dn ( dn \ ( 1 d \ dn ,. n . 

<M* 2 > 7TT— - • (42) 



d8-k L \d8-k L Jo \2P ks d5^ kL 1 2 7 <91ogo 8 
The point of this expression is that the non-gaussianity induces a bias, and that this can 
be expressed in terms of the unknown quantities (dn/d8- kL ) and dn/d logos that may be 
calculated in Newtonian gravity with (matter+A) only and /nl = 0. In the squeezed limit, 
the correction to dn/d5 kL from radiation approaches the value 8/35 ■ a eq log a eq sin 2 6 times 
dn/d logo's f° r P{k) oc k~ 3 in the approximation of ffl8l) . where 6 is the angle between k^ 
and k 2 . For isosceles triangles in the squeezed limit, 6 = 7r/2; the more precise derivation 
in [3H] demonstrates that 9 should be averaged over in the bias with measure sin#<i#. When 
ki > lOfceq, the contribution to the same quantity from /nl is 2f^/M kL times dn/dlogas, 
which is a significantly smaller contribution for large k and /nl = 1- However, the contribution 
from /nl is /c-dependent whereas the leading radiation contribution is not. Thus, the leading 
radiation contribution appears as just an overall rescaling of the galaxy bias and does not 
contaminate the signal for / NL when / NL is fit to match the shape of the galaxy power 
spectrum. Such contamination occurs only at sub-leading order and is comparable to /nl ~ 1- 
This can be seen in Figure where we show 5Q, which is proportional to the bias, after 
subtracting off the limiting value of 5Q at large k, and we compare it with the analogous effect 
from /nl = 1- Since the effect is order one, as we explained before, it becomes comparable to 
the GR effects that we expect to affect the distribution of collapsed objects at this order. 

5 Discussion 

Our main goal was to quantify the contribution from radiation to nongaussianities in the dark 
matter distribution compared with those from primordial nongaussianities. We have found 
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Figure 9: The corrections \5Q\ — |5<5fc>fe eq | for k\ = k<i = 10k 3 , compared to the contribution 
from primordial nongaussinities with /nl = 1. The correction is shown (solid, black) for the 
exact numeric result with fl mt o = 0.27, O n o = 3 x 10.~ 4 Q mt0 . The contribution from /nl is 
shown in blue, dashed. 



that the effect on the three-point function for modes with k « 10 k CCL is comparable to /nl ~ 4. 
At shorter wavelengths, the contribution from promordial nongaussianities shrinks like 1 / log k 
while the contribution from radiation becomes scale-independent. Nongaussianities may also 
be measured through the additional scale-dependence they induce on the power spectrum. 
Such scale dependence arising from radiation corrections is comparable to /nl = 1- We 
have also compared the size of primordial nongaussianities with the relativistic corrections in 
Newtonian gauge to the nonlinear growth of nongaussianities. Such corrections are negligible 
at large k ^> k cq but at smaller than k < k eq are comparable to /ml ~ 4 on equilateral 
triangles. However, at small k the effect on 5p/ p of a gauge transformation from Newtonian to 
synchronous gauge becomes large and parametrically the same as the relativistic corrections. 
Thus, to clearly interpret this result, Sp/p must be related to present-day observables. 

We have made a number of approximations to simplify the analysis. We neglected baryons 
and neutrinos, as well as higher order photon moments. We also have focussed on the dark 
matter density perturbations, which must eventually be related to visible objects. Further- 
more, the density field p(x) we have been using is the local inertial density, which will receive 
volume and redshift distortions from the local metric and peculiar velocity, as well as distor- 
tions along the line-of-sight. If observations of nongaussianities are to be interpreted with 
uncertainty less than A/nl < 4, then such distortions most likely must be understood and 
quantified as well. 
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Appendix 

A Second Order Sp/p in a Matter-Only Universe 

The form of the second order perturbations in a matter-only universe has been previously 
derived in [30j [31] using the formalism we adopt in this paper, and also independently in 
[34] using an action approach. The result can be written analytically in this case. For the 
convenience of the reader, we present and comment on the results in our formalism. 

In addition to setting Q r = 0, Q m = 1 in our equations of motion and ignoring the 
radiation perturbations, one must also choose appropriate initial conditions. We separate 
out the propagation from the initial conditions by using the method of Green's functions, as 
follows. First, it is useful to supplement the equations of motion we have so far with the 
component of Einstein's equation 

Hk 2 ¥ 2) + k 2 ^ - ^H 2 ikV® = s B (k, V ) = (43) 

2{k.q 2 )&^ + m 2 i{k.q l )V^5<£ . 

For reference, we present all components of the Einstein tensor for a metric in Newtonian 
gauge with only scalar modes: 

a 2 G° = -e~ 2 *3(ft - $') 2 + e 2 * ((<9$) 2 - 2«9 2 $) , (44) 
a 2 G^. = 8) [-e~ 2 * {{H - ®)(H - 3$' - 2tf') + 2(H' - $")) + e 2 * ((<9^) 2 + d 2 ^ - <9 2 $)] 
+e 2 ® [-diW?!! + diQdj® - 2<%$6> i) * - ^^(-$ + *)] , (45) 
G i0 = 2{H - $')<9^ + 2d& . (46) 

It is now straightforward to obtain a second order equation for 5^ sourced by the first order 
perturbations. We first use equation (I25p to eliminate = S4 + & 2 \ Then from (12411431) . 
we have 

, 2) 2k 2 S 3 -6HS 5 + 3k 2 H 2 5W +2k 4 &V 
V = 1 ' (47) 

Eliminating $'( 2 ) from equations f)26|27l) . we obtain $^ 2 ^ in terms of 5^ 2 \ 5 ,<y2 \ Si, which upon 
substitution back into (l26]) gives the desired equation for 5^ 2 \ Using that in a matter-only 
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universe, 7i = 2/rj, we obtain 

(2) „ 1296 + 72k V - 2fcV (2)/ 6 fc 2 (-42 + fcV) (2) 

216?7 + 18A; 2 r7 3 - A; 4 r? 5 216 + 18P^ 2 - fc 4 r/ 4 [ 1 



(48 



where .FfSy is a function of the source terms above, and therefore quadratic in the first order 
solutions. It is fairly long[f|, but we will see that the final answer simplifies considerably when 
we substitute the first order solutions. The above equation is a second order sourced equation 
of motion for 8^ and thus its solution is given by 

5(2) = c (?) 5l (r)) + c®5 2 ( V )+ P T[Si}{rf)G{ri,rf)drf , (50) 

Jo 

where G(rj,rj') is the Green's function 

[v,v) <W)<W) - <W)W) ' 1 } 

and Si, S 2 are the growing and decaying mode homogeneous solutions 

5,(7]) = 12 + (k V )\ (52) 
-18 + (fcry) 2 

*>(?,) = (M 5 ■ ( 53 ) 

The homogeneous solutions for 5^ are the same as the homogeneous (growing and decaying) 
mode solutions for and indeed for S^ 1 ' at any order, since their equations of motion only 
differ by the source term T . 

At first order, the equations of motion may easily be solved, and one finds that $ is 
constant and the matter perturbation growing mode is 

4 1) = -(12 + ,V)(|) , (54) 
Vt 1] = -(2ik V ) (|) . (55) 
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Explicitly, 

-l 



F[Si] = (fc 2 ry(fcV + 18) 2 (fcV-18fcV-216))" 



x 



(-rfS^k 12 + r? 9 S 3 (r7)fc 10 + rfS'Mk 10 + 1188t7 5 S^k 8 - 6t] 8 S 5 (t])k s 
-6r] S S' 3 (f])k s - 108ry 6 S*^(r/)fc 8 - 3f] 9 S' 3 \r])k 8 - 108it] 6 S' 2 {t])k 7 - 1188ry 5 5 3 (r/)fc 6 
+388877 3 S , 4 (?7)fc 6 - 540? ? 5 > S';(?7)fc 6 + 18r; 8 ^'(?7)fc 6 - 3888^ %{r])k i - 11664077S , 4 (77)fc 4 
+7128T] 4 S 5 {r 1 )k 4 ~ 3888t] 3 S[(r])k 4 + 9720r 1 4 S' 3 (r 1 )k 4 + 5832077 2 5 4 (ry)fc 4 + 1620t] 5 S 3 (t])k 4 
-it] (fc 8 ry 8 - 1188fcV - 3888fcV + 116640) S 2 (v)k 3 + 58320ir] 2 S 2 {r))k 3 
-2 (-fcV - 18/cV + 1620fc 4 77 4 + 25272/cV + 69984) Si(t])k 2 + ll6640r]S 3 (r))k 2 
+62208r) 2 S 5 (r])k 2 + V&Urr 2 S' 3 {rj)k 2 + 419904^(r;)fc 2 - SSSSOt? 3 ^)*; 2 - +166At] 3 S 3 {r])k 2 
-9720?7 4 S'5 (r])k 2 + 419904i,S 2 (77) k ~ 279936S , 5 (^) - 4199047755(77) - 69984t7 2 S , 5 ' (77)) (49) 
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Upon substituting these into the source term J 7 , the Green's function integral may be per- 
formed analytically: 



2\2 



+ 



23k 2 



V 

2 „2\2' 



167 ( 2 2 v 45 (? 2 - g 2 



14 



36 



where we have factored out $ 2 /36 in order to more easily compare the leading term with the 
PT kernel. To obtain the correct initial conditions, we take the early-time super-horizon limit 
r] — > 0, k — > as before (see paragraph preceding equation fl331 ): 

2 



<j(2) 



& - 25 4 - 2$ (2) 



3-^2 

^(2) = $ (2) + ^ 

Outside the horizon during the matter era, ( = — $ — = — 1$ — ^64, and thus $^ 2 ^ 



initially Using the limiting values of S3 
we obtain the initial condition for 5^: 



10 

A- 2 



(57) 
-§64 



(2) 



5fc 4_ 3(g 2_ g 2 ) 2 + 2A; 2 (g 2 + g 2 ) 

k 4 



91 92 - 



(5? 



Since the initial conditions are set at 77 — > 0, there is no component of decaying mode in 
equation (1501) . and thus ci = 5^ t /12. Putting everything together, we finally obtain 



1^(2) 

2°k 



9 (9 2 - 9 2 2 f 



59A; 2 
lT 



125 

u 



(9 2 + S2 2 ) 



+ (90 + 36 9l ^ 2 - 54 
V k 



7k 2 

fa 2 - & 

k 4 



4 1} (gl)^(g2 

36 



(59) 



(56) 



This agrees with the expression for 5p/p one may obtain from the second-order metric in [34]. 
We may rewrite this in a slightly different form to more easily compare with the PT result of 
equation (jSJ), which corresponds to the krj — > 00 limit: 



±5 (2) 



- a*) + • ft) 



91 92 

92 9i 



+ "fc(9i • 92) 2 + Ik 



X 



5, 



(1) 

9i 



5, 



(1) 
92 



1 + 37^1+3^ 

2 59ft 2 45ft 4 
7 + T4P + ~2F~ ' 



1 - 



ft 2 54ft 4 



2 A; 2 



fc 4 



7fc 



91 _ 92 

92 92 



m 2 9H 4 

+ 



2fc 2 



2fc 4 



(60) 
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A.l Squeezed Limit of 

We note one immediate check of equation (j591) based on the "squeezed" limit q\ ~ k ^> g 2 , at 
leading order. In this limit, is a very long wavelength mode whose only physical effect is 
to modify the physical time and wavelength of the short wavelength mode [H [36j [37]. More 
precisely, coordinates (77, x) with a long wavelength background $^ are equivalent to the 
modified coordinates (r]',x'), where 

a(r]')dx = (1 - $ kl )a(ri)dx , (61) 
a{rf)drf = (1 + $ h )a(v)dv ■ (62) 

During the matter-era, a{rj) oc r/ 2 , and thus 

7/ = 7/(1 + ^0 , (63) 
x' = x{l - jj$ fc; ) . (64) 

The effect of this coordinate redefinition on the short wavelength matter density p ks is 

1 dp ka 5 dp ka 
Pk a -> 0k,+ o$k,~j + o $ V 



3 l d\ogr] 3 '<91ogfc s 

1 5 OS 

Ph s + p-$ fei {-&iH6 k . + rj6' ks ) + -<$> h pk s -^ 



= p\ [ 5 + 2(72-6kW)[^^jj . (65) 

In terms of the explicit expression for 5^ (eq. [59]) . the effect of the long wavelength mode 
is obtained from — > (?i = 92 = + |^ (?i = ^,?2 — k s ) in the limit 
ki <^ k s , and we find 



i^- 2(72 -6kW) , (66) 



which indeed agrees with equation (jbol) . 
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